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In the hydrodynamic approximation we obtain analytic solutions to the Gross-Pitaevskii equation 
with positive scattering length which describe expansions of the Bose-Einstein condensates in quasi- 
one and quasi-two dimensional geometries. The expansion laws are expressed in terms of the initial 
sizes of the condensate and the trap frequencies before release, that is in terms of experimentally 
measurable parameters only. Three-dimensional effects are estimated with the use of variational 
approach. The analytical formulae show good agreement with available experimental data. 

PACS numbers: 03.75.Kk 



I. INTRODUCTION 

The properties of Bose-Einstein condensates (BECs) in lower-dimensional geometries, when the particles motion 
is tightly confined to zero point oscillations in one or two directions, have attracted much attention 'l'|-'l4'| because 
of specific coherence properties of degenerate Bose gases in lower dimensions and possible experimental realization of 
quantum gases beyond mean field approach. In order to realize a one-dimensional (ID) condensate, optical dipole traps 
have been used to achieve necessary asymmetry between the trapping frequencies (see, e.g. 0). Two-dimensional 
(2D) condensates have been created in an array of pancake-shaped traps formed by the periodic potential of a ID 
optical lattice 0- For large lattice potential depths, the motion along the direction of the array is frozen and the 
condensate is effectively split up into several smaller independent condensates confined in the individual lattice wells. 

Much experimental information on BECs in three-dimensional confining traps has been provided by imaging the 
density of the atomic cloud after switching off the trap potential. In the mean-field approximation, evolution of BECs 
in diluted gases is described by the time-dependent Gross-Pitaevskii (CP) equation jlj] 

Bih ft 2 

ih^ = - — /^i^ + V,^t{v)iP + g\i^\^, (1) 

where 

K.t(r) - \m{ujlx^ + + ulz^) (2) 

is the external trap potential, g = Anfi^as/m is the nonlinear interaction constant being the s-wave scattering 
length of the interacting atoms, and the condensate wave function ip is normalized to the number of atoms 

iV'pdr = N. (3) 

If the number of atoms is large enough, the CP equation can be reduced to the hydrodynamic equations, which admit 
simple self-similar solutions describin g as oscillations of BEC in the parabolic trap potential, so its free 3D expansion 
after switching off the trap potential |l5j|-[T8j. This theory shows excellent agreement with the experiment. 

In the recent experiments the free expansion of BEC has been investigated in ID waveguides [3| and in the array of 
2D pancake BECs 0| , and this line of investigations looks very promising. In expansion of BEC in a waveguide 
has been considered for quasi- ID case but without taking into account the quantum pressure in transverse direction. 
In I23 the effects of the quantum pressure have been calculated for stationary states with account of two transverse 
modes and expansion dynamic has been studied numerically. In stationary states of BEC in cigar-shaped and 
pancake traps have been considered by means of proper variational approach without discussion of free expansion 
dynamics. The aim of this article is to give analytic treatment of the lower dimensional expansion of BEC. To this 
end, under conditions that 3D CP equation can be reduced to its lower dimensional counterparts, in Section II we 
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solve these reduced equations in their hydrodynamic approximations with initial conditions corresponding to the 
equilibrium states of a confined in the parabolic trap BEC with positive scattering length (a^ > 0). Free expansion of 
BEC takes place either in longitudinal direction of the ID waveguide, or in transverse direction of the 2D pancake trap. 
If the conditions of lower dimensional reductions of GP equation arc not fulfilled, then effects of 3D hydrodynamic 
flow can be estimated with the use of variational approach what is done in Section III. Section IV is devoted to 
comparison of the theory with experiment. 



II. FREE ID AND 2D EXPANSION OF BEC 

The GP equations admits formulation in the form of the Hamilton least action principle with the action functional 

J Ldt, L = J Cdr, (4) 

where the Lagrangian density is given by 

£ = |(^:^ - VtV*) + I^IV^I' + VeM' + l9m\ (5) 

so that the corresponding Lagrange equation coincides with the GP equation (QJ. In cigar-shaped or pancake con- 
figurations it is not difficult to determine conditions when the tightly confined degrees of freedom become irrelevant 
("frozen") and the GP equation reduces to the effective ID or 2D equation. This problem has been considered in 
many papers (see, e.g. Chapter 17, and references therein) and we shall present here only the main points with 
the aim to define the relevant parameters and to formulate the applicability conditions. 



A. One-dimensional expansion 



If the axial frequency lu^ of the axially symmetric trap is much less than its transverse frequency lu± , 

A = ^ « 1, (6) 

UJJ_ 

and characteristic energy of zero-point transverse motion is much greater than the nonlinear energy, then the transverse 
motion reduces to the ground state of quantum oscillations in the radial trap potential. Characteristic radius of this 
motion is equal to a± = (fi/mwj,)^/^, and if we denote by Zq the characteristic axial size of the BEC, then with 
account of the estimate N ~ jV'pa^^o (see Eq. ©) the corresponding condition can be written in the form (see, e.g. 

mi) 

(7) 

If this condition is fulfilled, then we can factorize the condensate wave function 

i;{r,t) = ^{x,y)^{z,t), (8) 

where 

'/'(a;, y) = -1^ — cxp (-^^-Trr-) ' (9) 

is the ground state wave function of the transverse motion. Substitution of into Q,® and subsequent 

integration over transverse cross section lead to the expression for the action in terms of effective ID Lagrangian 
density 



CrD = f (*r* - *.**) + + \mu^lz^m'' + ^1*1' (10) 



(irrelevant constant term has been omitted). Now the evolution of \E'(z,t) is governed by Lagrange equation which 
yields the ID GP equation 

1 

ih-^t = "2;^*-^ + 2™^^^^* + (11) 



3 



where 



giD = ^ = ^^^"^ (12) 



is effective ID nonlinear coupling constant and ^' is normalized according to 

\^dz ^ N. (13) 



Equation Hll() governs ID longitudinal dynamics of the BEC in cigar-shaped trap. 
By means of well-known substitution 

*(z,<) = ^p{z,t)exp J v{z',t)dz'^ (14) 

equation Hll() transforms to the system 

Pt + (pv), = 0, (15) 

, 5lD , 2 , ( pI Pzz\ „ 

-*+-^ + — P^+-.- + ^(^i^-^j^-0- (16) 

The last term in Eq. called "quantum pressure" can be neglected, if it is much less than the nonlinear one 
proportional to .gizj, that is the condition 



is fulfilled. Then Eq. lfTC|) reduces to 



^ « — (17) 



Vt+VV:, + ^P:,+L0lz^Q (18) 

m 



which together with Eq. H15() provides the hydrodynamic approximation to description of BEC dynamics in parabolic 
trap potential. 

The stationary solution of the hydrodynamic equations gives the well-known Thomas-Fermi (TP) distribution of 
density for ID BEC: 

, , iN f z^\ 

where the integration constant has a meaning of the axial half-length of the condensate and can be expressed in 
terms of the number of atoms N , 

Zq = (3Ara,aiA-2)i/3, (20) 

Conditions (Q and p7|l of applicability of ID hydrodynamic approximation can be written upon substitution of l|2UI) 
as (see [ni) 

VA< ^<T- 21 

a± A 

Now we suppose that the axial trap is turned off, so that BEC expands freely in the axial direction. At the same 
time it remains being confined in the transverse direction, that is the radial motion is still "frozen" to the quantum 
ground state Hence, this expansion is governed by the hydrodynamic equations l|15() and 118() with the initial 
conditions H19(l . Similar problems were studied in nonlinear optics »22] for the opposite sign of the pressure term 
in (|18|l and for description of free 3D expansion of BEC jTEL ll^ FitL Il8l|. Here we use the same approach to the ID 
expansion of BEC in a "waveguide" . We look for the solution of Eqs. (|15|) . H18|) in the form 

P^^''^ = WoW){'-Zm)' -(-'^)--^(*)' (22) 
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where bz{t) and az{t) must satisfy the initial conditions 

6,(0) = 1, a,(0) = 0. (23) 
Substitution of (j221l into (^SIjCHI) yields the relation 

azit) = k{t)/hz{t) (24) 

and the equation for hz{t) 

k^u^l/hl (25) 

which can be integrated at once to give (see [2^, where this solution was applied to the description of the initial 
quasi-lD stage of expansion of pancake BEC in axial direction) 

V2L0zt = ^bz{bz - 1) + i In [26, - 1 + 2^6,(6, - 1)] . (26) 

This formula determines implicitly 6, as a function of time t. The function a{t) can be easily expressed in terms of 
bz{t), so that velocity field is given by 



The edge point of the density distribution moves according to the law 

Zmaxit) = Z(,bz{t) (28) 

with the maximal velocity 



Vma. (t) = = V2ZoLJz . / 1 - -i- . (29) 

at y 6, [t) 

In the limit of large t ^ uj^ ^ we find 

6,(t)~V2w,t, t>w7\ (30) 
so that distributions of density and velocity fields get simple form 

'iUmax t \ "max'' / 

and maximal velocity becomes constant 

t^maa; - V2ZqUJz, i > CJ^^ (32) 

These formulae describe a hydrodynamic flow "by inertia" when the density becomes so small that the nonlinear 
pressure does not accelerate the gas anymore. The formula (|32f) is convenient for comparison with experiment since 
the asymptotic value of the maximal velocity is expressed in terms of two experimentally measurable parameters — the 
axial frequency of the trap initially confining the BEC and the initial half- width Zq of the axial TF distribution. 
From H31|) we find the asymptotic distribution function of velocities 

, , , 3N / \ , , 

p[v)dv = 1 \dv, \v\<Vmax- (33) 

^Vmax \ ""max / 

The mean kinetic energy ( "release" energy) is given by 

m f 1 1 
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B. Two-dimensional expansion 



We can achieve 2D dynamics of BEC, if the axial frequency of the trap is much greater than its radial frequency 
Lu±, that is the inequality (jSJ is replaced by the inverse one, 



A = — > 1. 



(35) 



In this case we assume that the axial motion along the z axis is "frozen" , that is the energy of its zero point oscillatory 
motion with the quantum amplitude az — {h/muizY^^ is much greater than the nonlinear energy. This condition with 
the use of an estimate TV ~ \ip\'^RQaz, where Rq denotes a characteristic radius of the density distribution in the 
pancake (x, y) plane, leads to the inequality 



Na^^fRo 

az \ CLz 



If this condition is fulfilled, then we can again factorize the condensate wave function: 
where 

1 



{■A 



1/4 1/2 



exp 



z 
2^ 



(36) 



(37) 



(38) 



is the ground state wave function of the axial motion. Substitution of (|37|l . (|38|) into Q, 10 and integration over 
longitudinal coordinate lead to the expression for the action in terms of the effective 2D Lagrangian density 



ih , 

y 



re 

2m ' 



1 



9 



^ 1\J I'Kaz 



The corresponding Lagrange equation yields the 2D GP equation 

imt 

where = 9^ + 9^ is the transverse Laplace operator, g2B is the effective nonlinear interaction constant, 

a 



1 . . . 



92D = 



27ra, 



mo. 



and 5* is normalized according to 



\^^dxdy = N. 



(39) 



(40) 



(41) 



(42) 



The NLS equation (|4()|l governs 2D transverse dynamics of the BEC in a pancake trap. 
By means of a substitution 



*(ri,t) = Vp(rx,^)exp(^ / v(rl,t)-dr' 



(43) 



where = (x, y) and velocity v = {vx-, Vy) has components only in the transverse plane, the equation H4U|) transforms 
to the system 



Pt + Vi(pv) =0, 



(44) 



92D, 



vt + (vVi)v + — V_LP + wlr 



m 



2m? 



V 



4p2 2p 



= 0, 



(45) 
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where Vj_ = {dx, dy) is the transverse gradient operator. 

The quantum pressure term can be neglected if it is much less than the nonlinear term what leads to the condition 

1 < -. (46) 

Then Eq. (|45|l takes the form of the hydrodynamic equation 

Vt + (vV_L)v+^ViP + tjir = 0. (47) 
m 

The stationary solution of hydrodynamic equations corresponds to 2D TF approximation 

v^O, (48) 
where = + and the radius Rq of the density distribution is determined by the number of atoms N: 
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1/4 



i?o = \^-^Na,ai\'j . (49) 

It permits one to transform the inequality ()36|l to more convenient form, so that the applicability conditions of 2D 
approximation in the TF limit become 

1 < < A^. (50) 

When the radial trap is turned off, the BEC expands in the radial direction, but remains confined in the axial one. 
This radial expansion is described by the hydrodynamic equations H44|l and H47|l with the initial condition l)48|l. Now 
we look for the solution in the form 



.(r,i)=ra^(0, (51) 



where v is the radial component of velocity, and b± (t) and a±_ (t) must satisfy the initial conditions 

6^(0)^1, ai(0) = 0. (52) 

Substitution of (|51l) into (|44|) . H47|) leads to the relation a± = b±/b± between b±{t) and a±{t) and to the differential 
equation for 6_l(<): 

b^ = ul/bl. (53) 
Its integration with the initial conditions &_l(0) = 1, fc±(0) = 6^(0)q;j^(0) = gives 

b^{t) = ^l + Lolt\ (54) 

and hence aj^(i) — b±/b± — uj'^t/(l + uj'^t'^). Thus we obtain simple formulae for radial distributions of density and 
velocity fields: 

The edge point of the radial density distribution moves according to the law 



r„,a4t)=Ro^Jl+u}lt^ (56) 

with maximal velocity 

Vmax[t) = —7— = , ,) ^„ ■ (57) 

at yn + ujj¥ 
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In the limit of asymptotically large time t ^ uj^'^ we find 

2N 1 ^ \ r 

K'-'^)-^^^ vir,t)^- (58) 

I' Umax 1- \ Umax'' / '' 

where 

Vmax^Ro^^, i > tjj\ (59) 

As analogous formulae in ID case, these formulae describe a flow "by inertia" . Again the maximal velocity is expressed 
in terms of two experimentally measurable parameters — the initial radial frequency uj±_ of the trap before its switching 
off and the initial radius of the TF distribution. 
The asymptotic distribution of velocities is given by 

4A^ / v'^ \ 

p{v)dv = 1 - 3— vdv, (60) 



V \ V 

7nax \ max 



and the release energy is equal to 



III. EFFECTS OF 3D FLOW: VARIATIONAL APPROACH 



When conditions Q or H36|) are not fulfilled, the motion along the shortest size of the condensate cannot be 
considered as frozen one and we have to include into consideration the radial (for cigar-shaped trap) or axial (for 
pancake trap) flow. As was shown in 24j, it can be done by means of simple variational approach. However, in our 
case of strongly asymmetrical traps it is more natural for large enough number of atoms to assume the TF rather 
than Gaussian distribution along longer direction of the condensate, as it was done for static equilibrium case in |l3j| . 
Here we apply similar approach to dynamics of BEG. 



A. Cigar-shaped trap 



For cigar-shaped trap with TF distribution of density along the axial z direction we take the variational trial 
function in the form 



ip = A exp 




1 ^ exp 



* / 2 2 , 2 2\ 



(62) 



where w±, w^, a±, are assumed to be functions of time. It is supposed that Nas/Z^ ^ 1, that is in the radial 
direction the TF distribution is not achieved yet and therefore the radial wavefunction can be approximated well 
enough by the Gaussian function. The variable A is connected with widths w±, Wz by the normalization condition 
(O which yields 



A = 



1/2 



(63) 



Substitution of (|62|l into , ((SJ yields after integration the averaged Lagrangian 

fi^ da I 



L 

N ' 
daz 



2 dt 

at 

2m ^ 









1 2" 


2m 









1 



2 \2 



2 SNasti^ 1 



(64) 



5m 



W^Wz 



where we have neglected the term compared with since according to the estimates Wz ^ Zq, az ^ 

(m/h)wz/wz mujz/fi- ~ -^/^l the condition l/w^ ^ cez coincides exactly with the condition A^/^ ^ Nas/a±^ 
of applicability of the TF approximation in the axial direction (see (|21ll ). 
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The Lagrangian 1)64(1 leads to well-known formulas 



m 1 dw± m 1 dwz 

n w± at n at 



and to equations of motion for widths: 

ft2 GNa.h'^ 1 



w± + LJj^w± = 3 H — ^ 3 (66) 



2 2 3iVa,;i^ 1 

+ A ti^^Wz = 5 2 9 • ("') 

TO^ W\Wz 

These equations differ from those in 0] by numerical factors due to use of the TF distribution in the axial direction 
and by absence of the "quantum pressure" term in Eq. (|67|) which inclusion would be illegal in our approximation. 
In dimensionless variables 

fc± = — , bz = ^, T = LO±t (68) 



1 " ^ r7 







Eqs. (|Snil,(E3 take the form 



d'^b± , 1 2(7 1 



where the parameter 



(71) 



measures the strength of the nonlinear pressure in the radial direction. If g <C 1, then we can neglect the second term 
in the right-hand side of 169|) and obtain the stationary solution b±Q = b^o — 1 which corresponds to ID approximation 
of Section II. It is easy to see that in this case Eq. ((7())l reduces, for the case of free axial expansion, to ((25)1 . 
The equilibrium values of b± and bz are determined by the equations 

OzQ O OzO 

which differ from similar equations of |23| only through the notation. They determine the state of BEC before 
expansion. 

Expansion of BEC after switching off the axial trap is described by the equations 

dT^ ^ b\^ b b\hz' dr^ blbl' ^ ^ 

which can be easily solved numerically with the initial conditions 

b±iO) = b±o, j'±(0) = 0, 6,(0) = 6,0, &.(0) = 0, (74) 

where 6j_o and 6,o are determined by Eqs. 1(72(1 . Dependence of b± ~ w^/a^ and z^ax — Z^bzit) — Zbz{t)/bzo, where 
Z — ZobzQ is the initial half-length of the condensate, is shown in Fig. 1 for different values of the number of atoms 
N. For \/\a±/as <C <C a±/{Xas) we return to the analytical solution of Section 2. For TV > a±/{Xas) evolution 
starts from radial width greater than the quantum oscillator length aj_ and 6^ — )■ 1 for t — )■ oo. The axial expansion 
for T ^ 1 tends to the motion by inertia with constant velocities of atoms. It is easy to find the value of maximal 
velocity with the use of conservation law of Eqs. 1)73)) : 
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b |(t) 



z (t) [mm] 

max ^ ^ L J 




FIG. 1: Dependence of the dimensionless radial size b± — w±/a± of the condensate (a) and its axial size Zmax (b) on 
dimensionless time variable r — tujz during ID expansion in a cigar-shaped trap with the following values of the parameters: 
ax = 5 /im, as — 5nm, A — 0.05. The curves correspond to the numbers of atoms A'^ — 10^ (1), 10^ (2), 10^ (3). 



For r = we have the initial values (|74|) and for r — s- cx) we have dbz/dr — > bz^max, b± 

bz,max — '^^ 

Then with the help of Eqs. lO we find 



1, bz 



oo, hence 



b'iobzQ 



- - 7— 

q V o±o 



dz„ 



dt 



2Zujzbj_o 



(76) 



(77) 



For q 1, when b±Q = I, Z = Zq, we return to the ID formula H32|l . and for g 1, when b±o ^ 1, we obtain 
Vmax — ^UzZ. Thus, with increase of the number of atoms not only longitudinal size Z increases, but also the 
coefficient of proportionality between v,nax and uJzR increases from to 2. 

Asymptotic distribution of velocities as well as the expression for the release energy hold their form H33(I . H34|I . where 
Vmax is given by (|77|) . Dependence of the release energy on the number of atoms N is shown in Fig. 2. 



B. Pancake trap 



For pancake trap with TF distribution of density in the transverse direction, we take the variational trial function 
in the form 



■ exp 



2wV 



exp 



f 2 2 , 2 2\ 



(78) 



where the time-dependent parameters A, w± , Wz , a± , are connected due to the normalization condition ^ by the 
relation 



A 



2N 1 

7^3/2 uu'^Wz 



1/2 



(79) 



Substitution of (|78|l into , ^ yields after integration the averaged Lagrangian 



L 

N 



h? da± 
~2~dt 



duz 

2 dt 2m V w 



1 — i 




^ 2 \ 2 









w2 
2 



1 „ 

2m ^ 2 3 
8 Na^h^ 1 
3\/27r m Wz 



(80) 
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FIG. 2: The release energy in ID expansion of the condensate in a cigar-shaped trap as a function of the number of atoms A''. 
The energy is measured in units of the axial trap oscillator quantum huiz. 



where we have neglected the term l/w\ which is much less than due to the condition Nag/az S> 1 of applicability 
of the TF approximation in the radial direction. 

The Lagrangian H80|l leads to the relations H65|) and to equations of motion for widths: 

2 16 Nasfi^ 1 

W±+UJ^Wj_^—= 5 3 , (81) 



'2tt to^ wl^Wz 



..,2 2 16 Nastf 1 

Wz+X UJ±Wz = ^-^ + ^= 2 2-^, (82) 

m'^wi J2-K rn'^ w\ wi 



which in dimensionless variables 



o_L = — , bz = — , T^U!j_t, (83) 



where = aj_/\/X, take the form 



+ 6^ = 7T7-, (84) 



dT^ b\bz 



(Pbz , A2 A^q 1 



The parameter 



(86) 



measures the strength of the nonlinear pressure in the axial direction. If g ^ 1, then we can neglect the second term 
in the right-hand side of H85|l and obtain the stationary solution b^^ — bz — ^ which corresponds to 2D approximation 
of Section II. It is easy to see that in this case Eq. (|84|l reduces for the case of free axial expansion to (|53|l . 



11 



R (t) [mm] 
max L J 



b , (t) 




FIG. 3: Dependence of the radial size Rmax of the condensate (a) and its dimensionless axial size bz = Wz/az (b) on dimensionless 
time variable r = tuji_ during 2D expansion in a pancake trap with the following values of the parameters: Oz = 5 fj.m, Os = 5 nm, 
A = 20. The curves correspond to the numbers of atoms TV = 10^ (1), lO"* (2), 10^ (3). 



The equilibrium values of b± and are determined by the equations 

1 . Q 1 



-4 
LO' 



Ll6 



1 + ^-. 



(87) 



which again coincide with equations of '2l'| up to the notation. They determine the state of BEC before expansion. 
Expansion of BEC after switching off the radial trap is described by the equations 



1 



bib. 



6? 



VqJ_ 
"3" 616!' 



which can be solved numerically with the initial conditions (|74|) . where 6^0 smd b^o si's determined by eqs. (|87|l . The 
dependence of Rmax — Rb±{t)/b±Q and b^ — Wz/uz, where R = i?o&±o is the initial radius of the condensate, is shown 
in Fig. 3 for different values of the number of atoms N. For Qz/ag N azX^/ag we return to the analytical 
solution of Section 2. For N > a^A^/a^ the evolution starts from axial width Wzo greater than the quantum oscillator 
length az (i.e. bzo > 1) and 62 — *■ 1 for r — > 00. The radial expansion for r ^ 1 tends to the motion by inertia with 
constant velocities of atoms. The maximal radial velocity can be found with the use of conservation law of Eqs. (|88|l . 



dbj_ 
"d7 



which yields 



3 



dbz 
dT 



A2 



1 

6! 



b\bz 



— const. 



(89) 



b^^bzo ' 2q 



— I bzQ 



JzO 



Then for maximal velocity we obtain 



dR 



dt 



Ri 



I 1 + 36^0 



(90) 



(91) 



For q <C 1, when 620 — I, R = Rq, we return to the 2D formula H59|l . and for q ^ 1, when bzo ^ 1, we obtain 
Vmax — ■\/3/2ti>^i?. Asymptotic distribution of velocities as well as the expression for the release energy hold their 
form (|60|) . Ht)l| ). where Vmax is given by (|91|l . Dependence of the release energy on the number of atoms N is shown in 
Fig. 4. 
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2.10"^ 4.10"^ 6.10"^ 8.10"^ 10^ 

FIG. 4: The release energy in 2D expansion of the condensate in a pancake trap as a function of the number of atoms A^. The 
energy is measured in units of the radial trap oscillator quantum hij± . 

IV. DISCUSSION 

Here we compare the theory with available experimental data. 

Expansion of the condensate in quasi-lD waveguide was studied experimentally in Q]. In this experiment, the 
half-length of the condensate axial size was Z ~ 100 /im, and the time of transition to expansion by inertia was about 
20 ms which gives lOz ~ 50 s-^ Then formula ^ gives the value of the maximal velocity 

Vmax — V2ZuJz — 7mm/s 

which, in view of our approximate estimates of relevant parameters, agrees very well with the experimental value 
Vmax — 5.9mm/s. 

Free expansion of BEC in an array of 2D traps has been investigated in It is clearly seen in Fig. 2 of this 
paper that no axial expansion is observed what means that the condensate is effectively split into several independent 
condensates confined in the individual lattice wells. Therefore we can apply 2D theory to description of free radial 
expansion of each separate condensate. From Fig. 2 of ^ we find that maximal radial velocity is equal to v„iax — 
1.5 — 1.7mm/s, the initial radius is i? ~ 13/im, and the radial trap frequency is uj± ~ 132 s^^. Then formula l|59|) 
gives 

Vmax — -Rw_L — 1.7mm/s 

in excellent agreement with the above experimental value. Thus, as we see, the theory agrees with experimental 
data. Experimentally observed deviations from theoretical predictions could serve as indications on realization of 
non-mcan-ficld regimes of BEC behavior. 

In conclusion, we have considered in detail the theory of one- and two-dimensional expansion of BEC. Three- 
dimensional effects have been estimated with the use of variational approach. Good agreement with the available 
experimental data is found. 
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